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Abstract 

We experimentally demonstrate that supersaturated solution of sodium 
acetate, commonly called 'hot ice', is a massively-parallel unconventional 
computer. In the hot ice computer data are represented by a spatial 
configuration of crystallization induction sites and physical obstacles im- 
mersed in the experimental container. Computation is implemented by 
propagation and interaction of growing crystals initiated at the data-sites. 
We discuss experimental prototypes of hot ice processors which compute 
planar Voronoi diagram, shortest collision-free paths and implement AND 
and OR logical gates. 

Keywords: crystallization, Voronoi diagram, shortest path, unconven- 
tional computing, logical gates, physics of computation 

1 Introduction 

The last decade has witnessed a rise in the number of unconventional com- 
puting devices based on principles and mechanisms of information processing 
in physical, chemical and biological systems [11 . Most results in unconven- 
tional computing are theoretical [141 115] whilst there is a handful of intriguing 
experimental prototypes. The experimental laboratory research focuses on im- 
plementation of computing schemes in novel materials which helps us to grasp 
the inner nature of physics of computation. The following prototypes, or their 
families, give us a good representation of novel computing substrates. 

Reaction-diffusion computers are implemented in a spatially extended chem- 
ical media. Data and results of computation are represented by spatial distri- 
butions of reagent concentrations. Computation is realized by propagating and 
interacting diffusive or excitation waves pTJ. Experimental laboratory prototypes 



of reaction-diffusion computers solve a wide range of geometrical problems, op- 
timization and logical computation £Q . 

Physarum computing is based on adaptive foraging behaviour of Plasmod- 
ium of Physarum polycephalum. In Physarum computers data are represented 
by the distribution of attracting (food) and repelling sources (light). Compu- 
tation is implemented by the Plasmodium which optimizes it feeding pattern 
under control of attracting/repelling forces. Results of the computation are 
represented by the morphological structure of the plasmodium [25]. 126] [Jj [13] . 

Extended analog computers represent data as a spatial configuration of 
point-wise sources of current applied to a conductive medium. Computing op- 
erations are determined by physical properties of the conductive medium, (e.g. 
foam or gel). Results are represented by the distribution of voltages across the 
medium [24] . 

Other less-known prototypes include analog computing schemes using mi- 
crowaves [20], plane tessellation in gas-discharge semi-conductor systems [TU] , 
learning networks of memristors [31 , logical circuits in liquid crystals [22 , and 
glow discharge systems computing the shortest path [27] . 

Most experimental prototypes of unconventional computers either require a 
tailored hardware interface (analog computers, liquid crystals) and specialized 
equipment (memristors, gas-discharge systems), or they may have intrinsic lim- 
itations on the speed of computation (reaction-diffusion chemical processors). 
Physarum computer is the simplest to build but the most difficult to control, due 
to the sensitivity and somewhat unpredictable behaviour of the living creature. 

We aim, therefore, to provide an example of a novel computing material 
which is cheap to build, requires minimal resources to operate, implements 
computational procedures relatively quickly and is capable of solving a wide 
range of computationally-hard tasks. We show that sodium acetate trihydrate 
(colloquially called 'hot ice' due to its resemblance to ice and its crystalline 
behaviour) perfectly fits our specification of an 'ideal DIY unconventional com- 
puter' because it solves a variety of tasks by traveling and interacting waves of 
crystallization in its supersaturated solution. 

The paper is structured as follows. Experimental methods are outlined in 
Sect. [2] Section [3] presents experimental results on computing the Voronoi di- 
agram. In Sect. ^] we show how to extract direction towards the site of crys- 
tallization induction. Computation of one-destination-many-sources paths by 
crystallization of sodium acetate, and further extraction of a collision-free short- 
est path, are discussed in Sect. [5] Experimental implementation of logical gates 
is presented in Sect. [6] We summarize the findings and discuss further develop- 
ments in Sect. [7J 

2 Methods 

We prepare a supersaturated clear solution of sodium acetate trihydrate 
CH3COONa-3H20, pour the hot solution into Petri dishes and cool the so- 
lution down to -5°C. To induce crystallization we briefly immerse aluminum 
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wire (powdered with fine crystals of the sodium acetate) into the solution. 

To compute shortest paths in a space with obstacles we mimic obstacles with 
blobs of silicon, small Petri dishes fixed to the bottom of experimental container 
and labyrinths made of Blu-Tack®[^] 

Dynamics of crystallization is recorded with FujiPix 6000 digital camera 
and still high-resolution images are produced by scanning containers with crys- 
tallized sodium acetate in HP Deskjet 5100 scanner. Magnified images of the 
crystalline structure formed are taken using Digital Blue® QX-5 USB micro- 
scope. 

Software tools for extracting directions of crystal growth from still images 
and models of cellular-automaton shortest path calculation are coded in Pro- 
cessinEpl 



3 Voronoi diagram 

Let P be a nonempty finite set of planar points. A planar Voronoi diagram 
of P is a partition of the plane into such regions, that for any element of P, 
a region corresponding to a unique point p contains all those points of the 
plane which are closer to p than to any other node of P. A unique region 
vor{p) = {z s R 2 : d{p,z) < m)Vm g R 2 , m ^ z} assigned to point p 
is called a Voronoi cell of the point p. The boundary of the Voronoi cell of a 
point p is built of segments of bisectors separating pairs of geographically closest 
points of the given planar set P. A union of all boundaries of the Voronoi cells 
determines the planar Voronoi diagram: VD(P) — U pe pdvor(p) [28| [23], 

A basic concept of parallel approximation of Voronoi diagram by wave pat- 
terns propagating in a spatially extended medium [5] is based on time-to- 
distance transformation. To compute a bisector separating two given points 
p and q we initiate wave-patterns in the medium's sites geographically corre- 
sponding to p and q. The waves travel the same distance from the sites of 
origination before colliding. The loci where the waves meet indicate sites of the 
computed bisector. 

Precipitating reaction-diffusion chemical media are proved to be an ideal 
computing substrate for approximation of the planar Voronoi diagram |30U18l [Tj. 
A Voronoi diagram can be approximated in a two-reagent medium. One reagent 
a is saturated in the substrate, drops of another reagent (3 are applied to the 
sites corresponding to planar points to be separated by bisectors. The reagent 
(3 diffuses in the substrate and reacts with reagent a. Colored precipitate is 
produced in the reaction between a and f3. When two or more waves of diffusing 
a meet, no precipitate is formed |18j . Thus uncolored loci of the reaction- 
diffusion medium represent bisectors of the computed diagram. 

Most reaction-diffusion chemical processors are quite slow in computing a 
Voronoi diagram at the macro-scale, due to speed limitations caused by diffu- 
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sion. Crystallization patterns in sodium acetate propagate much quicker and 
thus could improve speed limitations of existing reaction-diffusion algorithms. 

We represented data points of set P by a configuration of pins (pieces of 
aluminium wire) fixed through the lid of a 9 cm Petri dish (Fig. [T|l) . The pins 
were powdered with fine crystals of sodium acetate. 

To start computation we place the lid on the dish with the supersaturated 
solution. The pins become immersed into the solution. They induce crystal- 
lization. Patterns of crystallization propagate - as classical target waves - from 
the sites of crystallization induction. 

A crystallization wave stops when it encounters another wave of crystal- 
lization (Fig. [ljac). Boundaries between stationary patterns of crystallization 
represent edges of planar Voronoi diagram VD(P). The experimental bound- 
aries perfectly correspond to bisectors of the diagram computed by the classical 
Fortune's sweeplinc algorithm (21] (Fig.JTJi). 

4 Detecting directions towards sites of crystal- 
lization induction 

Crystals growing from from nucleation sites bear distinctive elongated shapes 
expanding towards their proximal ends (Fig. |2^). Not only a crystal's overall 
shape but also the orientation of saw-tooth edges indicate the direction of the 
crystal's growth (Fig. [2J)). 

The direction of crystal growth can be detected by conventional image pro- 
cessing techniques, e.g. edge detection procedures (Fig. [2), or by a complemen- 
tary method of detecting directional uniformity of image domains as discussed 
below. 

We detect local directions of crystal grows in the following way. A fine 
grid of nodes is applied onto given image and for every node x of the grid 
we calculate a set V of vectors of length r originating at x and directed at 
angles a € {0,0.1,0.2, • • • , 2/tt}. For each vector v € V we calculate a sum of 
standard deviations a(v) of RGB colors of the image pixels coinciding with v. 
A vector v' = min^gv c( v ) with minimal sum of color deviations indicates a 
local direction of crystallization propagation. A configuration of local direction 
vectors, r = 20, is shown in Fig. [2ji. 

The vector configuration represents sinks determined by points from set P. 
If we place a mobile agent at any site but bisectors of the vector configuration 
Fig. [2}i the agent will be attracted to the closest data point p. This leads us to 
another problem solvable with sodium acetate — computation of a collision-free 
shortest path in a space with obstacles. 

5 Computation of a collision-free path 

Previously [3] we have designed a cellular-automaton algorithm for constructing 
one-destination-many-sources shortest path. The automaton is a regular array 



4 




(c) (d) 



Figure 1: Computation of Voronoi diagram by crystallization of sodium ac- 
etate: (a) configuration of pins representing data points P, (b) photograph and 
(c) scanned image (sites where crystallization was induced are encircled) of crys- 
tal structure induced by the pins, (d) planar Voronoi diagram of P computed 
by classical algorithm. See video of the experiment at [H] . 
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Figure 2: Detecting directions towards sites of crystallization induction: crystal 
pattern in xlO (a) and x60 (b) magnification, (c) edge detection operation is 
applied to image Fig. [TJi, (d) sample configuration of local direction vectors 
computed from image Fig. [TJ;. 
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Figure 3: Computation of shortest paths to a single destination: (a) example 
of pointer orientation (resting cells are gray, excited cell is marked with '+' 
and refractory cells with '-'), (b) example of pointer configuration in presence 
of obstacles and three shortest paths (bullets) to the site of initial excitation; 
orientation of cell pointers is shown explicitly. 



of finite state machines (cells) which takes discrete states and update their 
states simultaneously and in discrete time depending on the states of its closest 
neighbors. Here we consider hexagonal array of cells. Each cell has six neighbors 
and takes three states — resting, excited and refractory. A resting cell becomes 
excited if it has at least one excited neighbour. An excited cell switches to a 
refractory state, and refractory cell returns to resting state independently on 
states of its neighbors. When we excite one cell of a resting array a wave of 
excitation is initiated. The excitation wave propagates in the array. 

To store the computed path we supply every cell with a pointer. Pointers in 
all cells are set to nil initially. When a cell becomes excited its pointer orients 
towards the direction from where the excitation wave came. An example is 
shown in Fig. [3^. 
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Obstacles are imitated by always-resting, or non-excitable, cells. Exam- 
ples of the pointer configuration and shortest collision-free paths are shown in 
Fig. ^p. A cell assigned to be a destination (in the north-west corner of the 
array Fig. [3Jd) is excited. Waves of excitation propagate from the site of initial 
stimulation. The waves spread around obstacles and disappear in the absorbing 
boundaries of the array. The excitation waves induce orientation of cell pointers 
(Fig. ^jp) . Array of pointers represents a set of shortest collision-free paths from 
any cell (source) of the array to the site of initial excitation (destination) . Three 
such paths are shown in Fig. [SJd. 

As we demonstrated in Sect. [4] longitudinal crystals growing in sodium ac- 
etate are physical analogs of pointers indicating from where — locally - the 
wave of crystallization came from. That is crystallization of sodium acetate is 
the physical implementation of the original cellular-automaton algorithm j3] for 
computation of shortest collision-free path. 

Let us consider experimental realization. We represent obstacles by impen- 
etrable barriers and induce crystallization at a single site of the medium. For 
demonstration purposes we have tested various types of objects to represent 
obstacles, including 35 mm Petri dish glued with silicon to a bottom of larger 
Petri dish (Fig. |4^) and blob of silicon (Fig. [5J)) to make simple obstacles, and 
Blu-Tack® to make more complex obstacles (Fig.jZj;) and a labyrinth (Fig, [ijl). 

From the images (Fig. |4]) we calculate a configuration of local vectors (Fig.|5|. 
A vector at each point indicates the direction from where the wave of crystalliza- 
tion came from. The vectors represent an attracting field, which 'pulls' virtual 
objects, placed in the field, towards the destination (site of initial disturbance). 
We found that the attracting field is approximated by the vectors satisfactory 
for a wide range of vector parameters: vector grid steps 5 (Fig. ISb), 10 (Fig. |5fc,) 
and 20 (Fig.^d), and vector length 20 (Fig.[|i), 20 (Fig.[5]ib) and 50 (Fig.|5f). 

To check whether vector configurations extracted from scanned images of 
crystallized sodium acetate provide sufficient representation of many-sources- 
one-destination shortest paths we employed pixelbots [6]. A pixelbot is a virtual 
pixel-sized automaton. The pixelbot moves along the direction of nearest non- 
zero vector. The pixelbot imitates a tiny (2-5 mm) in diameter robot which 
travels along sodium acetate crystals or in the longitudinal grooves between the 
crystals. If there are no non-zero vectors in the pixelbot 's neighborhood the 
pixelbot moves at random [6]. 

Examples of test-run trajectories of pixelbots traveling in the vector con- 
figurations are superimposed on scanned images of crystallized sodium acetate 
in Fig. [6] Pixelbots easily reach the destination in the case of simple obstacles 



(Fig. |6pb) . They spent more time wandering at random in arenas with complex 

obstacles (Fig. [6J;) . Pixelbots do not perform well in a maze (Fig.[6]i): in some 
parts of the maze pixelbots navigate satisfactory and travel for substantial dis- 
tances (see trajectories in the western half of image Fig.[6]i), in other parts they 
become stuck forever falling into indefinite loops of repetitive motion (eastern 
half of image Fig. [6ji) . 
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(c) (d) 



Figure 4: Scanned images of crystallization of sodium acetate occurred in shal- 
low dishes with impenetrable obstacles: (a) 35 mm Petri dish glued to the bot- 
tom of the container, (b) blob of silicon, and complex obstacles (c) and maze (d) 
made of Blu- Tack®. Sites of crystallization induction (destinations of shortest 
paths) are encircled. See videos of the experiments at [8]. 



9 




10 




11 



6 Implementation of AND and OR logical gates 



We implement logical gates in a geometrically constrained, T-shaped, containers 
with sodium acetate. Shoulders of a T-shaped container represent inputs x and 
y and vertical stem — output z. When we induce crystallization in any one or 
both shoulders the whole solution in the T-shaped container, including z-output, 
becomes crystallized. Assuming liquid phase represents logical True and solid 
phase logical False, sodium acetate in T-shaped container implements OR gate. 

If we induce crystallization in both shoulders of the container crystalliza- 
tion waves collide at a junction and form a boundary in result of the collision 
(Fig. [7]). If crystallization is induced only in one of the shoulders no bound- 
ary is formed in z-output. This hints to the following interpretation. Pres- 
ence/absence of crystals in shoulders x and y represent True/False. Pres- 
ence/absence of boundary (bisector of crystallization induction sites in x and 
y) represents z =Truth/False. In such interpretation sodium acetate in T- 
shaped container implements and gate. 

7 Discussion 

We experimentally demonstrated that crystals growing in supersaturated so- 
lution of sodium acetate ('hot ice') compute planar Voronoi diagram, approxi- 
mate set of collision-free shortest paths in a space with obstacles and implement 
Boolean conjunction and disjunction in a geometrically constrained medium. 

Every micro- volume of the sodium acetate solution can input external stim- 
uli/data, participate in crystal growth and represent a boundary between col- 
lided waves of crystal growth. Data are inputted in the sodium acetate solution 
in parallel, e.g. we use a set of pins to input planar point set to compute a 
Voronoi diagram. Results of computation are also read in parallel, e.g. using 
optical recording. A micro-volume of the solution can be removed or added 
prior to computation, this is the sign of fault-tolerance. Therefore supersatu- 
rated solution of sodium acetate is a massively-parallel fault-tolerant computer. 

With regards to Voronoi diagram computation the hot ice computer is a 
highly accelerated version of precipitating reaction-diffusion chemical Voronoi 
processors [301 HE1 EZ1 HH US]- We demonstrate in experiments that bound- 
aries between crystal patterns perfectly represent edges of Voronoi diagram 
constructed on sites of crystallization induction. We did not show however how 
exactly the edges can be extracted from still images of the crystallized sodium 
acetate. This may be a topic of further studies. 

In laboratory experiments we proved that hot ice computer approximates a 
set of many-source-one-destination collision-free shortest paths in a space with 
obstacles. The shortest path problem is a typical benchmark task for uncon- 
ventional computing devices. There experimental evidences that shortest path 
can be computed in excitable chemical medium (Belousov-Zhabotinsky reac- 
tion) [5J [55] , in precipitating chemical medium interfaced with cellular automa- 
ton [5], and in plasmodium of Physarum polycephalum [2"5] , 
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(c) 

Figure 7: Experimental implementation of and gate in sodium acetate: (a) orig- 
inal scanned image of the crystallized sodium acetate, (b) contrast enhanced 
image, (c) scheme of the computation, two sites of crystallization induction 
(corresponding to True values of input variables x and y) are encircled, bisec- 
tor separating x and y is shown by a line. 
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These implementations, although very appealing, are however far from per- 
fect. Some of them require substantial assistance of external hardware/software 
in the detection of the 'intersection' of waves initiated in source and destination 
sites [HI [53] or extraction of a shortest path [5], others just do not produce a 
directed path at all [25 . Hot ice computer physically represents a many-source- 
one-destination set of shortest paths, where destination is a site of crystallization 
induction. If there would be a mobile robot compatible in size with width of 
crystals, such a robot would find the destination from any site of crystal pattern 
just by following the orientation of crystals. However, our experimental imple- 
mentation is not ideal: it is sometimes possible to take a wrong route or fall 
into cyclic route by following individual crystals. We will perfect the realization 
in future. 

We demonstrated that two logical gates - AND and OR - are implemented 
in crystallization of sodium acetate solution. This construction complements 
experimental realization of XOR gate in precipitating reaction-diffusion chemical 
medium [3] . One drawback of our approach is that different properties of the 
medium represent values of output Boolean variables: liquid and solid phases 
in OR gate, and presence/absence of a bisector in AND gate. Thus, in principle, 
AND and OR gates cannot be cascaded. One can consider three valued variables, 
where liquid phase represents 0, crystal phase - c and bisector (a boundary 
between collided waves of crystallization) - b. Then a T-shaped container with 
sodium acetate implements the following binary composition ★: 

* I b c 
c c 
b c c c 
c c c b 

Whether such operation is useful for implementation of logical circuits or not 
will be studied elsewhere. 

Finally we will say a few words about re-usability. In [1 we provide compar- 
ative analysis of existing non-linear medium computers in terms of re-usability. 

All excitable processors, including Belousov-Zhabotinsky medium, are re- 
usable chemical computers, because when development of excitation patterns is 
finished (and a dissipative structure of excitations representing result of com- 
putation recorded) the medium can be forced back to its resting state, thus 
being ready for next cycle of computation. However, excitable processors are 
memoryless, because they do not, as a rule, store results of computation in a 
non-volatile memory. 

Precipitating chemical processors, e.g. [3], have memory. They represent re- 
sults of computation in the spatial distribution of precipitate which, when dried, 
can stay unchanged for years. However, the precipitating chemical processors 
cannot be re-used. They are not re-usable but "disposable" . 

Hot ice computer is an experimental example of re-usable non-linear medium 
processor with non- volatile memory. Patterns of crystallization, which represent 
results of computation, stay almost intact for weeks and months. As soon as 
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they are heated above 58°C the patterns are transformed into liquid phase, thus 
being ready for next cycle of computation. 
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